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1 Introduction 



There have recently been several papers on hitting probabilities for systems of 
stochastic partial differential equations (SPDEs). The first seems to be [12], 
which studied mainly polarity of points for the Gaussian random field which is 
the solution of a system of linear heat equations in spatial dimension 1 driven 
by space-time white noise. Next, the paper [7] studied hitting probabilities 
for a nonlinear system of (reduced) wave equations in spatial dimension 1, 
and established upper and lower bounds on hitting probabilities in terms of 
Bessel-Ricsz capacity. 

The paper [3] considered a system of nonlinear heat equations in spa- 
tial dimension k — 1 with additive space-time white noise and established 
lower and upper bounds on the probability that the solution {u(t, x), {t, x) e 
M+ X [0, 1]) hits a set A C M'^ in terms of capacity and Hausdorff measure, 
respectively. These results were extended to systems of the same heat equa- 
tions but with multiphcative noise in [4]. The paper [5] extends these results 
to systems of non-linear heat equations in spatial dimensions k > 1, driven 
by spatially homogeneous noise that is white in time. Some other results on 
hitting probabilities for parabolic SPDEs with reflection are contained in the 



The objective of this paper is to begin a similar program for systems of 
stochastic wave equations, starting with the analogue of [3]. We note that 
properties of solutions of stochastic wave equations in spatial dimensions 
k > 1 are often much more difficult to obtain than their analogues for heat 
equations, due to the greater irrregularity of the fundamental solution of the 
wave equation. One example of this is the study in [8] of Holder continuity 
of sample paths for the 3-dimensional wave equation. 

In [3], various conditions on the density of the random vector 
{u{t, x), u{s, y)) were identified that imply upper and lower bounds on hitting 
probabilities. These were applied there essentially to study the linear case, 
and then the non-linear case with additive noise was handled by appealing to 
Girsanov's theorem. Because of the absence of a suitable Girsanov's theorem 
for heat or wave equations in spatial dimensions A; > 1 (a problem also noted 
in [5]), we will develop first some general results that will also be useful for 
nonlinear equations, and then apply them here to linear wave equations in 
spatial dimension k > 1, driven by spatially homogeneous noise that is white 
in time. More precisely, we consider the d-dimensional stochastic process 
U = {{ui{t, x),i = 1, . . . ,d),{t,x) e [0, T] X R'^} which solves the system of 



papers [20], [21], [6]. 



SPDEs 




d 



{t,x) - Aui{t,x) = ^(7ijF-'(i,x), 



(1) 
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for (t, x) g]0,T] X M^', with initial conditions 

UiiO,t)^-^{0,x)^0. (2) 

Here, A denotes the Laplacian on R*^, and a = (cTij) is a deterministic, 
invertible, d x d matrix. The noise process F := {F^, . . . , F^) is a centered 
(generahzed) Gaussian process whose covariance is informally given by an 
expression such as 

E{F\t,x)P{s,y) = dijd{t-s) \\x-y\\-^, (3) 

where 5ij denotes the Kronecker symbol, S{-) is the Dirac delta function 
at zero and (3 > 0. More precisely, let C(j^(]R^^^) denote the space of in- 
finitely differentiable functions with compact support, and consider a family 
of centered Gaussian random vectors F = {F{(p) = {F\(p), . . . , F^{(p)), (p E 
C^CUJ^'^^), with covariance function 

E{F{^)F{i;)) = / dr [ r{dx) {<p{t, •) * i;{t, ■)){x), (4) 

where ip{t,x) := ilj{t,—x) and F is a non-negative and non-negative defi- 
nite tempered measure on Ml'. We note that (4) reduces to (3) if V{dx) — 
\\x\\~^dx. By the Bochner-Schwartz theorem (see [16]), there exists a non- 
negative tempered measure on M*^ (termed the spectral measure of F) such 
that r = jT/i, where denotes the Fourier transform. Elementary properties 
of the Fourier transform show that equation (4) can be written 

E{F(^)F{i;)) ^ [ dr f a{di)J^^{t,.mmtrm- (5) 

JR+ JR'= 

Let G{t, x) be the fundamental solution of the wave equation. Generically, 
the solution u (1) is given by 

Ui{t, x) = J^^ G{t -r,x-y)Y, ^ijM^idr, dy) (6) 

where M = {M^, . . . , Af^) is the martingale measure derived from F (see [2] 
for details). However, it is well-known that G is a function in dimensions 
k e {1,2} only, so the stochastic integral in (6) should be interpreted in the 
sense of [1]. We note that according to (4) and (5), 

E{{u,{t,x)r)^(j:al] f^dr j^^^) |^G(t-r)(0|% 
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and it is well-known (see [18]) that 

^GWK) = ?^«. (7) 

Following [1], [13], we note that when ^ is not the null measure, the 
solution u{t,x) of (1) is a random vector, and the right-hand side of (6) is 
well-defined, if and only if the following hypothesis is satisfied: 

(H) 0</-i^<oo. 

In this case, the process u given by (6) is a natural example of an anisotropic 
Gaussian process as considered in [19]. Notice that for the covariance density 
in (3), condition (H) is satisfied when P G ]0, 2 A A;[. 

In Section 2 of this paper, we develop several results on hitting probabil- 
ities that are extensions of those of [3] and are appropriate for studying the 
wave equation in all spatial dimensions. Indeed, the results of this reference 
were taylored to the particularities of the heat equation in spatial dimension 
1, while our results highlight the role of the spatial dimension. Theorem 
2.1 gives a lower bound on hitting probabilities. Proposition 2.3 and Theo- 
rem 2.4 give upper bounds. These three results apply to arbitrary stochastic 
processes, while Theorem 2.6 gives a refinement of the upper bound in the 
case of Gaussian processes. These results are used in Section 4, but will also 
be useful for studying non-linear forms of (1), which is currently work in 
progress. 

In Section 3, we give simple conditions on a Gaussian process {X{t)) 
that ensure an upper bound on the density function of {X(t), X{s)). This 
is related to a result in [7]. The upper bound is expressed in terms of the 
canonical metric of the Gaussian process. 

In Section 4, the main effort is to obtain upper and lower bounds on the 
behavior of the canonical metric associated with the process u (Proposition 
4.1). This is somewhat intricate for the lower bounds, mainly because the 
expression for E{{u{t,x) —u{s,y)y) involves integrals of trigonometric func- 
tions, and these are not so easy to bound from below by positive quantities. 
Section 4 ends with applying the results of Sections 2 and 3 on hitting prob- 
abilities to obtain Theorems 4.4 and 4.5. These yield the following type of 
bounds: 

cCap._2(fc+i)(A) < P{u([to,T]x[-M, M]'')nA 0} < C n,_2(k+i) (A), (8) 

" 2-13 " 2-/3 

where Cap^ and denote capacity and Hausdorfi^ measure, respectively 
(their definitions are recalled in Section 2) . We note that the same dimensions 
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appear on both the left- and right-hand sides of (8). This conclusion could 
also have been deduced from Theorem 7.6 in [19], which contains a general 
result on hitting probabilities for anisotropic Gaussian processes. This is 
because our estimates on the canonical metric of u mentioned above, together 
with our Lemma 3.2, verify Conditions (CI) and (C2) in this reference. We 
note also that these estimates hint at the fact that condition (C3') of [19] 
should be satisfied by u. 

We recall that a point 2; e R'^ is polar for u if for all to > and M > 0, 



Notice, as a consequence of (8), that ii d < 2{k + l)/(2 — then points are 
not polar for -u, while ii d > 2{k + l)/(2 — f3), then points are polar for u. 
In the case where (3 is rational and 2{k + l)/{2 — (3) = is an integer, then 
polarity of points in the critical dimension d is an open problem. 

As mentioned above, in work in progress, we plan to extend these results 
to systems of nonlinear stochastic wave equations with additive noise but 
without using Girsanov's theorem. It is a separate endeavor to develop the 
estimates from Malliavin calculus needed for multiplicative noise, as was done 
in [4] for the heat equation, and which will also make use of the results in 
Section 2. 

2 General results on hitting probabilities 

Throughout this section, V = {v{x), x G M"^}, m G N*, denotes an M'^-valued 
stochastic process with continuous sample paths. We will fix a compact set 
/ C M"^ of positive Lebesgue measure and consider an arbitrary Borel set 
C R''. Our aim is to give sufficient conditions on the stochastic process V 
which lead to lower and upper bounds on the hitting probabilities 



in terms of the capacity and the Hausdorff measure of A, respectively, of a 
certain dimension. Here, v{I) denotes the image of / under the (random) 
map X I— >• v{x). 

We now introduce some notation and recall the definition of capacity and 
Hausdorff measure. For any 7 e R, we define the Bessel-Riesz kernels by 



P{z^u{[to,T]x[-M,Mf)} = Q. 



P{v{I) n A 7^ 0} 



if 7 > 0, 




(9) 



if 7 < 0, 
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where c is a constant whose value will be specified later in the proof of 
Lemma 2.2. Then, for every Borel set A C W', we define V{A) to be the set 
of probability measures on A. For fj, e V{A), we set 

K^{\\x-y\\)Mdx)ii{dy). 

J A J A 

The Bessel-Riesz capacity of a Borel set A C R*^ is defined as follows: 

Cap^(A) = 

with the convention that l/oo = 0. 

The 7-dimensional Hausdorff measure of a Borel set ^4 C K*^ is defined 
by H^{A) = oo if 7 < 0, and for 7 > 0, 

n^{A) = liminf |f;(2ri)^ : A C U~i5,,(a;,), supr^ < e\ . (11) 

Here and throughout the paper, Br{x) denotes the open Euclidean ball cen- 
tered at X and with radius r. Positive constants will be denoted most often 
by C or c, although their value may change from one line to the next. For a 
given subset S gMJ^ and u > 0, we denote by S^'''^ the i/-enlargement of 5". 

We begin by studying the lower bound for P{v{I) f] Ay^ $}. 



inf S^(u) 



(10) 



Theorem 2.1 Fix N > and assume that the stochastic process V satisfies 
the following two hypotheses: 

(1) For any x,y E I with x ^ y, the vector {v{x),v{y)) has a density Px,y, 

and there exist 7, o; e ]0, oo[ such that 

Px,y{zi,Z2) < Cj ^-^exp 

llx — ylp 

for any Zi, Z2 G [—N,NY, where C and c are positive constants inde- 
pendent of X and y. 

(2) One of the next two conditions holds: 

(P) The density px of v{x) is continuous, hounded, and Px{w) > for 
any x e I and we [-(A^ + 1), + l]'^. 

(P') For any compact set K gM.'^ and for any x e I, 'm.iyj^KPx{w) > 
Co > 0. 



C\ Zi - Z2\ 



\X 
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Then, for all Borel sets A C [— A^, N^, there exists a positive and finite 
constant c — c{N, a, 7, /, m) such that 

P{v{I) n ^ ^ 0} > c Cap^(^_^)(A). (12) 

Proof. Without loss of generality, we may assume that Cap 2.(^_^)(A) > 0, 
otherwise there is nothing to prove. Under this assumption, we necessarily 
have ^(7 — m) < d and A 7^ (see [11], Appendix C, Corollary 2.3.1, p. 
525). " 

Assume first that A is a compact set. Following the scheme of the proof 
of Theorem 2.1 in [3], we consider three different cases. 

Case 1: 7 — m < 0. Let z E A, e E ]0, 1[, and set 



{2s) 

We will prove that E{J^{z)) > ci and E[(Js{z)Y] < C2 for some positive 
constants Ci, 02- With this, by using the Palcy-Zygmund inequality ([11], 
Chapter 3, Lemma 1.4.1), and noticing that Cap^(^) = 1 for /? < 0, we 
obtain 

p{Mz) > 0} > > c 

= C Capi(^_^)(^). 

But P{Je{z) > 0} is bounded above by P{v{I) n A^'^ ^ 0}. Since A is 
compact and the trajectories of v are continuous, by letting e tend to 0, we 
obtain (12). 

The lower bound for E{Js{z)) is a direct consequence of assumption (2). 
To obtain the upper bound for E[{Js{z)y], we first use the hypothesis (1) to 
obtain 

EUizm <C [dx [ dy ^ 

Ji J I \\x — y\y 

Let po > be such that / C -Bpo (0). Fix x E I; after the change of variables 
y X — y and by considering polar coordinates, we easily get 



E[{j,,t{z))']<c r p'^-'-'^dp. 

Jo 



The last integral is bounded by a finite positive constant c(m, 7, /) Therefore 
we obtain E[{J^{z))^] < C2. 
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Case 2: < ^(7 — m) < d. Let 11 e V{A). Let = j^jd^Bs{o) and 

= I j^, /^(^^) 1^,(0) {v{x) - z) 

^ j^dx{ge* 1^) {v{x)) . 

Clearly, assumption (2) implies that E {J^{nj) > ci, for a constant ci which 
does not depend on /j, or e. Moreover, 



E[{Je{l^)y] = f dx [ dy [ dzi [ dz2 

Jl Jl JBe[z) JBe[z) 

X {ge * IJ) {Zi) {Qe * //) (^2)Px,j/(^l, Z2). 



By hypothesis (1), Lemma 2.2 below and Theorem B.l in [3], this is bounded 

by 

C , dzi dz2 {ge * IJ) {zi) {g^ * jj) {z2)K2, A\\zi - Z2\\) 

JBe{z) JBeiz) ' 

= C {ge * 

< C 

By choosing ii such that < 2/Cap2^^_^^(74), we obtain 

C 



E[{Je{fi)Y] < 



Cap2(^_^)(^)' 



and this yields (12) by a similar argument as in Case 1. 

Case 3: ^ — m — 0. The proof is done exactly in the same way as for 
Case 2, by applying Theorem B.2 in [3] instead of Theorem B.l. 

Now let A be a Borel set included in [—N, N]^. It is well-known that 

Cap^(A) = sup Cap^(F), (13) 

FcA, F compact 

(see for instance [9], Chapter 3). Therefore, for any compact set F C A, we 
have 

P{v{I) n A ^ 0} > P{v{I) n F 7^ 0} > c Cap2(,_„)(F). 

This yields (12) by taking the supremum over such F and using (13). 

The proof of the theorem is complete. □ 

In order to end the study of the lower bounds, wc prove a technical 
lemma which was used in the proof of Theorem 2.1 to relate joint densities 
with Bessel-Riesz kernels. 
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Lemma 2.2 Fix a, j e]0, oo[. There exists a constant C := C{N, a, j, I, m) 
such that for any a & ] — N,N[, 

Proof. Fix po > such that / C Bpsi{0). Fix x E I and consider the change 

of variables z — a~'^{x — y). Denoting by X the left-hand side of (14), we 
have 

X < cm a-^(^-™) / dz ^ exp ( | . 

Let 

J ^ dz -n— n- exp --n— n— ■ 

JB^iO) \\z\\^ V Ikll"/ 

Using polar coordinates, we have — J^i + with 

J, exp (-i), 

J.=/f exp (-1. 

Clearly J7i < C{pq,N). In order to study J'2, we bound the exponential by 
1 and we consider three different cases. 

Case 1: If m — 7 < 0, then 

^ <C(7V,a,7,Po,m). 

Case 2: If m — 7 > 0, then 

^ <C(7,Po,m)a^(^-). 

Case 3: If m — 7 = 0, then 

2, fN\ 

J2<-l0g - . 

Q! V a / 

Since X < C {I)a~i^'^~'^^ J' , wc reach the conclusion using the definition of 
Kj3{a) for /9, a G M (see (9); in the case where m — 7 = 0, the constant c in 
(9) must be chosen sufficiently large). □ 
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We now study upper bounds for the hitting probabihties. For this, we fix 
5 > 0, £ e ]0, 1[, ji, . . . , e Z, and we set j = (ji, . . . , j^) and 

m 

^i=nM'0^+iM]- (15) 

1=1 

The next statement is an extension to higher dimensions of Theorem 3.1 
in [3]. 

Proposition 2.3 Let D C and 7 > 0. We assume that there exists a 
positive constant c such that, for all small e e ]0, 1[, z e D^^\ and any set 
R^j such that R]f\I 

P[v{R])f\B,{z)^%}<ce'. (16) 

Then, there exists a positive constant C such that, for any Borel set Ad D, 

P{v{I)nA^%} KCn^-r^^A). (17) 

Proof We suppose 7 — y > 0, otherwise 7i-y_^(A) = 00 and therefore (17) 
obviously holds. Clearly, by the additive property of probability, 

P{v{i)r\B,{z)^%}< P{v{R'^)r\B,{z)^^], 

for any e > 0. Since / is bounded, the number of terms in the sum on the 
right-hand side of this inequality is bounded by a multiple of e^T. Hence 

P {v{I) n B,{z) ^ 0} < Ce-fp {v{R'.) n B,{z) ^ 0} . 
Then, using to (16), we obtain 

P {v{I) n B,{z) ^ H}} < Ce^-T . (18) 

This yields (17) by a covering argument, as it is shown in the proof of Theorem 
3.1 in [3]. For the sake of completeness we sketch this argument. 

Fix e G ]0, 1[ sufficiently small and consider a sequence of open balls 
{Bn, n > 1) with respective radii r„ e ]0,s], such that B^ n A ^ 0, A C. 
[Jn>iBn and 

Y,{2rnr-f<n,-:.{A)+e. 

n>l 
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Then, by (18), 

p {vii) nA^<D}<Y.P {v{i) n s„ ^ 0} 

n>l 
n>l 

< C {n,.^{A) + e) . 

Finally, we let £ | to conclude. □ 

In the next theorem, we give sufficient conditions on the process V for the 
assumptions of Proposition 2.3 to be satisfied and therefore to ensure (17). 

Theorem 2.4 Let D C W'. Assume that the stochastic process V satisfies 
the following two conditions: 

(1) For any x e , the random vector v{x) has a density p^, and 

sup sup Px{z) < C. 

(2) There exists 5 £ ]0, 1] and a constant C such that, for any q e [l,oo[; 

Ei\\vix)-v{y)r)<C\\x-yt. 

Then for any 7 G ]0, d[, the inequality (16) holds, and consequently, for every 
Borel set Ad D, 

P{v{I)nA^{b} <Cn^-ru{A). (19) 

Proof. We keep the notations of Proposition 2.3 and write = (jie^,l — 
1, . . . , m). For any z G D^^^ and Rj such that Rj D I 0, set 

Yf = \\v{x') - z\\, Z] = sup ||v(x) - v{xm. 

By applying the version of Kolmogorov's criterion as it is stated in [14], 
Theorem 2.1, page 26, using assumption (2), we obtain 

E{{z^y)<c\\x-x^^\r, 

for any q G [1, oo[ and a G ]0, — ^[. Hence 

E {{Z^Y) < C£^°^ (20) 
10 



with 7o < 1 - ^. 

Let 7 e ]0, d[. We first prove that 

P{Z] > \y^} < Ce\ (21) 
For this, we consider the decomposition 

and then give upper bounds for each term on the right-hand side. 

Clearly, from the boundedness of the density stated in assumption (1), 

P{r/ <£^} <C£T, 

and by Markov's inequality along with (20), 

■' 2 

Therefore, 

P{Z'j > Wf} < C [e< + , 

for any 7o < \ — Since 7 G ]0, we can choose 70 < 1 and q arbitrarily 
large such that ^ < 70 < 1 — ^. Hence, we obtain (21). 

If v{Rf} n B^{z) ^ 0, then < £ + Therefore, 

P{v{R'^ n B,{z) 7^ 0} < P{F/ < e + Z^} 

< P{Z| > iy/} + P{F/ < 2£} 
<c(e< + s*^) 

since 7 e ]0, This proves (16) for any 7 G ]0, d\. According to Proposition 
2.3, we obtain (19). □ 

The remainder of this section is devoted to extending the validity of (19) 
to 7 = (i in the case where V belongs to a particular class of Gaussian 
processes. For this class, we will prove that, instead of (20), the following 
stronger property holds: 
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For any e G ]0, 1[, for each j G with fl / 7^ 0, and every q G [1, cxd[, 
there is C > such that 



E sup \\v{x)-v{xm'i\ < C£^. 



(22) 



Then, we will show that P{v{Rfj n B^{z) 7^ 0} < Ce'^ (see Theorem 2.6 be- 
low). Together with Proposition 2.3, this will yield the desired improvement. 

We first give a sufficient condition which applies to arbitrary continuous 
stochastic processes V. 

Lemma 2.5 Let u G ]0, 1]. Suppose that for any e G ]0, 1[ sufficiently small, 

\v{y) -v{y)\\ ' 



E \ I dy dy 

lBe{x) JBeix) 



exp 



\y-y\ 



< Ce 



2m 



(23) 



where C is a positive constant. Let S^{x) — {y E R"^ : — y|| <£"}. 

Then, for any q G [1, oo[; there exists C > such that for all small £ > 0, 



E{ sup \\v{x) -v{y)\\'^] < Ce'^. 



(24) 



Proof By (23), B{uj) < 00 a.s., where 



B(u) 



dy / dy 

s^{x) Js-{x) 



exp 



\v{y){uj) - v{y){uj) 

\\y - y\V 



We apply the Garsia-Rodemich-Rumsey lemma (cf. [17], exercise 2.4.1, pg. 
60) to the functions ^^{x) = — p{x) = x^ and functions / : S^{x) C — > 
given by the sample paths of the process V restricted to the parameter 

set >S'^(a;), to obtain 



\v[x) — V 



-n^-yW ^ (CiB{oo)\ 



u 



where Ci is a positive constant which depends only on m. Consequently, for 
any q G [l,oo[, 

2.^ ^_,(C,B{uj)\ 

' ' vu du. 



E ( sup \\v{x) - v{y)\\A <SE\ r xjj 
\yes^{x) I \Jo 



u 



2m 



We notice that since ip ^{x) = In(l + x) is an increasing function on [0, 00), 
the constant Ci above can be taken arbitrarily large. In the sequel, we will 
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fix g G [1, cxd[ and we will assume that Ci > (e^ ^ — 1) C2 ^2'^"\ where C2 is 
the square of the volume of the unit ball in R"*. Then 



B{uj) > €28— > 



_ 1 



2m 



for any u e [0, 2£^]. 

Jensen's inequality applied first to the convex function fi{x) = x^, a; G M 
and the integral with respect to the measure ii{du) = u^~^du, and then to the 
concave function (p2{x) = ln''(l + a;), x G [e^~^ — 1, oo[ and to the expectation 
operator, yield 



1 

2£l7 



£;( sup \\v{x)-v{y)\A <%e'^~^ 

JO 

Jo 



E 



InM 1 + 



C^Biuj) 



u 



2m 



v!' ^du 



1 

2e7 



InM 1 + ^^^^ 1 vu'-^du, 



u 



2m 



with C3 = CiC. With the change of variable m — > — , we have 



2fi 



2m 



^0 V 



£ 

Ce. 



dw 



This proves (24). 



□ 



We can now sharpen the result of Theorem 2.4 in the case of Gaussian 
processes. 

Theorem 2.6 Assume that the stochastic process V = {v{x), x G M™} 
is continuous, Gaussian, centered, with independent, identically distributed 
components {vi{x), x G M™}, i = 1, . . . ,d, and inf^gj(i) Var(vi(a;)) > 0. Fix 
5 G ]0, 1] and suppose that for any e > small enough and any Rj (defined 
in (15)) such that R]f\I 



e(I dy I dy 



exp 



\v{y) -v{y)\\ 
\\y - y\\^ 



<Ce s . 



(25) 
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Then for every z eM.'^ and Rj as before, 

P{v{R'j) n B,{z) ^ 0} < Ce''. (26) 

Consequently, for any Borel set A C W^, 

P{v{I) n ^ 7^ 0} < CHd-^iA). (27) 

Proof By Lemma 2.5, assumption (25) implies (22). We use this property 
and adapt the proof of Proposition 4.4 of [3]. First, for any z e W^, we write 

P{v{m n BJz) 7^ 0} = P i inf \\vix) - z\\ < e\ . 

yxeR^. J 

Next, we write the condition \\v{x) — z\\ < e in terms of two independent 
random variables, as follows. Set 



E {vi{x)vi{xf)^ 



Because y is a Gaussian process, 

E {v{x)\v{x''jj^ = c^j{x)v{x^j). 

Set 

Yj = inf^ \\c'j{x)v{x'j) - z\\, Zj = sup \\v{x) - c^(x)w(a;^) 
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Again because y is a Gaussian process, these two random variables are in- 
dependent, and 

P{ inf \\v{x) - z\\ <e}< P{Y' <e + Z'}. (28) 

Our next aim is to prove that for any r > 0, 

P(r/ < r) < Cr^ (29) 
For this, we first notice that by independence of the components of V , 

P{y!<r)<i{P{Gl,), 

i=l 

where 

xeR 



= { M \c'^{x)v,{x'^) - Zi\ < r}. 
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By setting = infj-gj^e C^jix), we have 



Since V is centered and infj.gj(i) Var(i>i(x)) > by hypothesis, Schwarz's 
inequahty and (22) yield 



- 1| 



E 



Vl{Xj) (vi{x) - Vi{Xj)'j 



< C 



Var(fi(a;^)) 

E ( \vi{x) -vi{x'j) 
YariviixD) 



2\ \ 2 



(30) 



for any x e This imphes ^ < Cr, and since the density of Vi{x^j) is 
bounded, we get 



p(v,(x^j)eB^(z)^<Cr, 



and therefore (29) holds. 

By (29) and the independence of and ZJ, 



P{Y' <e + Zn < CE 



Consider the decomposition = Zj'^ + Zj'^, where 

Z."'' = sup ||^;(a;) - v(a;^)||, = ||^;(a;^)|| sup |1 - c'Ax)\. 



By (22), we have E 



< Ce'^. Moreover, by (30) and (25), 



E 



7£,2 



< 



Ce'^E {\\v{x'j)\\'^) < Ce'^. 



This completes the proof of (26). Finally, (27) follows from Proposition 2.3. 

□ 



15 



3 Joint densities of Gaussian processes 

Consider a Gaussian family of centered, R'^-valucd random vectors, indexed 
by a compact metric space (T, rf), that we denote by X = {Xt, t E T). 
We suppose that the component processes {X^, t G T), i = l,...,d, are 
independent. We also assume mean-square continuity, that is, by letting 

5{s,t)^{E{\\X,-X,r))'' 

denote the canonical (pseudo)-metric associated with X, we have S{s,t) 
as d{s, t) — > 0. 

Let Ps,t{^i: Z2) denote the joint density of {Xs, Xt) at {zi, Z2) G M^*^. The 
purpose of this section is to establish upper bounds of exponential type for 
4(^1,^2)- We notice that these conditions, and in particular, condition (c) 
below, are easily verified in many examples. 

Proposition 3.1 Suppose that 

(a) af- := Var{Xl) > 0, for any i = 1, . . . ,d and for all t eT, 

(h) Corr(X], XI) < 1, for any i — 1, . . . ,d, s,t eT with s ^t, 

(c) there exists rj > and a positive constant C > such that, for all 
s,teT, 

sup Wli-<^\<Ci^i'^t))'^'■ 

i&{l,...,d} 

Fix M > 0. Then, there exists C > such that for all s,t E T with s ^ t 
and zi,Z2 G [-M,M]'^, 

for some positive and finite constants C and c. 

Proof. Notice that (a), (b) and the independence of the components yield 
the existence of Ps,t- 

Fix i — 1, . . . , d, and denote by p\ ^{zl.| Z2), p^g{-\z2) and pl{-) the joint 
density of {XI, XI) at {zi,Z2), the conditional density of XI given XI = Z2 and 
the marginal density of XJ, respectively. It is well-known (linear regression) 
that 
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where 



2 2/1 2 \ 

^s,t^(^t{^- Ps,t), Ps,t 



and, for the sake of simpUcity, wc have omitted the index i. Since 

the triangle inequahty along with the elementary bound (a — 6)^ > — 6^ 
yields 



PL(^1,^2) < 



2'KUsTs-i 



exp 



U^l - ^2 



4r2 



X exp — 5 exp 



By hypotheses (a) and (c), s h- >• a1 is bounded above and bounded below by 
a positive constant, therefore for Z2 G [— M, M], 



P',,t(^l,^2) < exp 



exp 



'M2|1 - ms,t\ 



The conclusion now follows from Lemma 3.2 below and the independence of 
the components oi X. □ 

Lemma 3.2 With the same assumptions and notations as in Proposition 
3.1, there exist constants < ci < C2 < oo such that for all s, t e T, 

(1) Ci 5{s,t) < Ts,t < C2 5{s,t), 

(2) \l-msA < C2 5(s,t). 
Proof. A simple calculation gives 



^2^2 2 



(31) 



(see [12], equation (3.1)). Therefore, by hypothesis (c) of Proposition 3.1, 



1 - Pit < 



C 



s{s,ty 



s'-'t 



Prom assumptions (a) and (c), it follows that there is a positive constant 
C2 < oo such that for all s, t e T, 

tI<cI 5{s,t)\ 
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which proves the upper bound in assertion (1). 

For the lower bound in (1), we note that for s near t, the second factor on 
the right-hand side of (31) is bounded below since 6{s, — > as d{s, t) — > 0. 
Further, by hypotheses (a) and (c), 

>Cl5(s,^)^ 

for s near t. This proves the lower bound in (1) when 5{s,t) is sufficiently 
small. 

In order to extend this inequality to all s, i e T, it suffices to observe that 
by hypothesis (b), 

if s 7^ t, and by hypothesis (c), for e > 0, there is c' > such that cr^ag—ag^ > 
c' for S{s,t) > e. This proves the lower bound in assertion (1). 
In order to prove assertion (2), observe that 



|1 - ms,t\ 

and 

< 5{s,tf + 5{s,t)at < c 5{s,t). 
This completes the proof. □ 



4 Hitting probabilities for the stochastic 
wave equation: the Gaussian case 

In this section, we consider the solution to equation (1), which is the d- 
dimensional Gaussian random field defined by 

u{t, x) = [ G{t -r,x- y)aM{ds, dy), {t, x) e [0, T] x (32) 

Since a is invertible, we may assume as in [3] that a is the identity matrix. 
Notice that, in this case, u{t, x) — {ui{t, x), . . . , Ud{t, x)), with 

Ui{t,x)=f*f G{t -r,x -y)M'{ds,dy), (t, a;) G [0, T] x M\ 
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1, . . . ,d, and therefore, the component processes {ui{t, x), {t, x) G [0, T] x 
, i — 1, . . . ,d, are i.i.d. 

Most of the results of this section require the following hypothesis: 

(H^) The spectral covariance measure /i is absolutely continuous with re- 
spect to Lebesgue measure on IR'^ and its density is given by 

/(0 = lieir'+^ Pe]0,2Ak[. 

Equivalently, r{dx) = C{k, I3)\\x\\~'^dx (see [10]). Notice that (H^) implies 
(H). 

In the sequel we fix a strictly positive real number to- We first aim for 
lower bounds on hitting probabilities. For this, we intend to apply Theorem 
2.1. The required upper bound on the joint densities will be obtained by 
combining Proposition 3.1 and the next two results. 

Proposition 4.1 Assume (Hp). Fix M > 0. Then, there exist positive 
constants Ci, C2 such that, for any {t,x), {s,y) e [to,T] x [—M,M]'', 



Ci{\t-s\ + \\x-y\\f-'' <E(\\ut,,-u 



<C\{\t-s\ + \\x-y\\f-r (33) 

Proof. The structure of this proof is similar to that of Lemma 4.2 in [3], but 
the methods for obtaining the estimates differ substantially. Without loss 
of generality, we will assume that d — 1. Let R{x) — E{u{t, x)u{t,0)) with 
t>to. Then 

E {{ut,, - ut,yf) = 2 (i?(0) - R{x - y)) . 

Following the steps of the proof of Remark 5.2 in [8], with the dimension 
/c = 3 replaced by an arbitrary value of k, and therefore the Riesz kernel 
||^||-(3-/3) replaced by we obtain 

i?(0) - R{x) < C\\xf-^. (34) 

We next fix y e R'^ and consider increments in time. Let to < s < t < T. 
Using (5) and (7), 

E [{u{t, y) - u{s, y)f) = S^{s, t) + ^2(5, t), 

with 

di |sin((t-r)||e||)-sin((s-r)||e||)P 



Sr{s,t)= r dr f 

Jo JM 

S2{s,t)^ f dr [ 

Js JM. 



* r di sm\{t-rm\) 
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With the changes of variables r — > s — r and C ^ {t — s)^, along with the 
trigonometric formula sin x — sin y = 2 sin cos , we obtain 

<C\t-s\^-f'. 

For the term 6*2(5, t), wc consider the changes of variables r ^ t — r and then 
^ — > r^, which easily yield 

<C|f-5|'-''. 

Hence, we have proved 

E ({u{t, y) - u{s, y)f) < C\t - ^l^-^^, (35) 

with a positive constant C depending only on T. With (34) and (35), we 
have established the upper bound in (33). 

We now prove the lower bound in (33) using several steps. 

Step 1. Assume s = t > to and x ^ y. The arguments in the proof of 
Theorem 5.1(a) in [8] can be trivially extended to any spatial dimension k. 
Therefore, there is a positive constant Ci such that for any x,y & [—M, M]^ , 

E({u{t,x)-u{t,y)f) >c,\x-y\^-^. (36) 

Step 2. We show that, for arbitrary x,y e [-M, M]*^ and to < s <t <T, 
E {iu{t, x) - u{s, y)f) > c\t - s\^-^. (37) 
Indeed the left-hand side of this inequality is equal to 
Ri{s,t;x,y) + R2{s,t; x,y), 

with 

R,{s,t;x,y)^[drlj^^ 

X \J^G{t -r,x- -m - ^G{s -r,y- , 
R2(s,t;x,y)^ dr J^^,^^\J^G(t - r,x - . 
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Since R2{s,t] x,y) is positive, we can neglect its contribution. (We notice 
that 

R2{s,t;x,y) > C|i - 

for some positive constant C. For k — 3, this is shown in the proof of 
Theorem 5.1 in [8], and it is easy to check that the arguments go through 
any dimension.) 

By developing the integrand in Ri{s,t]x,y), wc find 

Ur - r, X - -m - ^G{s -T,y- ■m\' 

= sin((t-r)||e||)-e^^-(^--)sin((.-r)||e||)f 
_l-cos(2(t-r)||^||)^l-cos(2(.-r)||^||) 



2 2 
- cos(e • {y - x)) [cos((t - sm\) - cos{{t + s - 2r)||e||)] . 

After integrating this last expression with respect to the variable r, we obtain 
a positive quantity which is the sum of the following three terms: 

Ai^s[l- cos{{t - s)m) cos(e • {y - x))] , 

= (cos(e ■ {y - X)) - cos((t - s)\\m , 

sin(2(t-.)||e||) Mit-sMW) 

As- 411^11 2|ieii '''^^ 

For the integration with respect to the variable ^, we restrict the domain to 
the set 

Do = {^eR': UUt -s)>l, cos((t - sMW) > 0}. 
Note that on Dq- wc have Ai > 0. In fact, 

A = 5 [1 - cos((t - s)U\\) + cos((i - - cos(e • {y - xM 

>s[l-cos{{t-s)m)]- 
Moreover, 

2 

1^2 + ^31 < 



Thus, with the change of variables i ^ {t — s)^, we easily obtain 



-slt-sf-f' [ „ ri-cos(||^| 

Alk||>l; cos(||«,||)>0| W r"^+^ 



dw 

l{\\w\\>l; cos(||«)||)>0} \\w\ 
>C2\t-s\'-^. 
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Similarly, 



All 



:\A2 + As\ < C3\t-S 



Therefore, by the triangular inequality, we obtain 



- q|3-/? 



Ri{s,t;x,y) > C2\t - s 



^-^-C3|t-.|3-^>||t 



|2-/3 



if \t — s\ < This proves (37) for small values of |t — s|. 

To extend the validity of (37) to arbitrary values of |i — s|, we notice 

that Ri{s,t; x,y) is a continuous and positive function of its arguments and 
therefore, it is bounded below on {(s, t; x, y) G [to, T]^ x [— M, M]^'' : \t — s\ > 
s} by some constant c^, for any s > 0. Hence, if 2T > |t — s| > we also 
have 

i?i(s, t; X, y) > c\t — s\'^'~^, 
for some c sufficiently small. 

Step 3. Suppose |t — <s| > 



Cl 2-/3 



4C2. 



x — y\, where ci appears in (36) and 



C2 in the right hand-side of (33). By Step 2, we clearly have 
E {{u{t, x) - u{s, y)f) > c\t - 

> c 



t — S 1 / Cl 

+ 



2-/3 



2-/3 



^2 '2 V4C2, 
>Cs{\t-s\ + \x-y\f-f'. 



x-y\ 



Step 4- Suppose |t — s| < 



Cl 



4C2 



2-0 



X — y\. Then, 



E ({u{t, x) - u{s, y)f) > ]-E ({u{t, x) - u{t, y)f) - E ({u{t, y) - u{s, y)f) 



2 

>\cAx-yr^-C2\t-sr^ 



> — 

- 4 



Cl (\x- y\ 



+ 2 



4C^ 

Cl 



2-/3 



\t-s\ 



>C,{\t-s\ + \x-y\f-^ . 
With this, the lower bound in (33) is proved. 



□ 
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Remark 4.1 (a) As mentioned in the Introduction, Proposition 4.1, to- 
gether with Lemma 3.2, estabhshes conditions (CI) and (C2) of [19] for the 
process U. 

(b) A consequence of the preceding proposition is that the sample paths 
of (32) are Holder continuous, jointly in {t,x), of exponent 7 G ]0, ^^[, but 
they are not Holder continuous of exponent 7 > ^y^. We refer the reader to 
[8] for a similar result on the solution to a nonhnear stochastic wave equation 
in spatial dimension k — 3. 



The next proposition is a further step towards proving that the process U 
satisfies the assumptions of Proposition 3.1. We denote by cr|,j. the common 
variance of Ui(t, x), i = 1, . . . ,d. 

Proposition 4.2 Assume that condition (H) is satisfied. Fix {t,x), {s,y) e 
[to,T] X R''. Then, 

(i) al,>C{toAtl)>0, 

(ii) \<jl^-(Tly\ < C\t-s\. 

If, in addition, we assume that for k' < k, all k' -dimensional submanifolds 



o/R are sets with null ^-measure, then 
(Hi) for any {t, x) 7^ (s, y) and i = 1, . . . ,d, 

Corr (u\t,x),u\s,yj^ < 1. 

Proof. The variance of u{t, x) is 

2 /* / sin^((t-r)||e||) 

= Jo L ^^^^^ — Mr — ' ^ ^ 

and satisfies 

C{t A t') < al <C{t + t') (39) 

(see for instance [15], Lemma 8.6). This proves (i). 
Assumption (H) implies that 

W(r|15||) 



SUP LM<)^™^<a 



re[0,T]JR>' Ur 

Hence, assuming to < s < t < T , we obtain 



iieip 
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which yields the conclusion (ii) of the proposition. 

We now prove (iii) by checking that for any (t, x) ^ (s, y) in [to, T\ x R'^, 

2 2 2 r\ 

where (Tt,x;s,y denotes the covariance of Ui{t,x) and Ui{s,y) for any i — 
l,...,d. 

Case 1: s < t. If crlx'^ly — <^tx-sy "were equal to zero, then the random 
variables u\t,x) and u'^{s,y) would have correlation equal to 1; therefore, 
there would be A e R such that u^{t,x) — \u^{s,y) a.s., and, in particular, 
we would have 

E[{u\t,x)-\u\s,y)f^^Q. 
The left-hand side of this equahty is 

fdr [ n{dO \J'G{t-r,x--mf 

Js JR* 

+ /' dr [ u(dO \J^G(s -r,x- - ^^G(s -r,y- , 

Jo JR'^ 

which is bounded below, as in (39), by C((t — s) A (t — s)^). This leads to a 
contradiction. 

Case 2: s = t, x y. We start as in the preceding case by assuming that 
^Ix^^ly - ^t,x;t,y = and hence 

E(^(u%x)-Xu%y)f^^O 

for some A e R. The left-hand side is equal to 

r dr I jx{di) e^«- - Ae'«-^ ' |^G(r, ^(Of ■ 

Jo jRk 

If A = 1, then the integrand vanishes when cos[^ • (x — y)] — 1 or sin(r||^||) = 
0, which occurs on a (/c — l)-dimcnsional manifold of R'^. Hence, by the 
assumption on fi, we reach a contradiction. 
If A 7^ 1, then 

f'dr [ a{dO e'^-^-Xe'^-y'\J^G(r,-mf 

Jo Jr'' 



> fdr [ a{dO (1-A)2 



sin^(rllCll) 



This last integrand vanishes only when sin(r||^||) = 0. Thus we also get a 
contradiction in this case. The proof of the proposition is now complete. □ 

We can now obtain the required properties on densities, as follows. 
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Proposition 4.3 Assume (H^). Fix M,N > and {t,x), {s,y) e [to,T] x 
[-M,M]*^ with {t,x) 7^ {s,y). 

(a) Let pt,x;s,y{- , ■) denote the joint density of the random vector 
{u{t,x),u{s,y)). Then 

Pt,x;s,y{^l,Z2) < 3(5^^ exp -— — , 

{\t - s\ + \x - y\)^ V {\t - s\ + \x - y\) 

(40) 

for any Zi,Z2 G [—N, N]''', with C and c positive constants not depending 
on (t,x), (s,y). 

(b) Let pt^x denote the density of the random vector u{t,x). Then, for each 

{t, x) e [to, T] X M'^ and z e [-N, N]'^, 

PtM > C, (41) 

and 

sup sup Pt,x{z) < C. (42) 

ze[-N,N]'i {t,x)e[to,T]xR'' 

Proof By Propositions 4.2 and 4.1, we see that the process U satisfies the 
hypotheses of Proposition 3.1 with rj — Thus, we have statement (a). 
The density pt^x is given by 



2 ) ' 

t.x I 



with as in (38). By (39), we obtain both (41) and (42). □ 

The next theorem gives lower bounds on hitting probabilities. 

Theorem 4.4 Assume (H^). Let I, J be compact subsets of[to,T] and R*^, 
respectively, each with positive Lebesgue measure. Fix N > 0. Then, we have 
the following: 

(1) There exists a positive constant c — c{I, J, N, (3, k, d) such that, for any 
Borel set A C [-N, N]"^, 

P{u(I X J)nA^$}>c Cap^_ 2jk+ix (A) . (43) 

" 2-/3 

(2) For any t E I , there exists a positive constant c — c{J, N, /3, k, d, t) such 

that, for any Borel set A C [-A^, A^]"', 

P{ii(0} X J) n 0} > c Cap^_^(A). (44) 
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(3) For any x E J, there exists a positive constant c = c{I , N, P, k,d,x) 
such that, for any Borel set A C [—N, N]''', 

P{u (/ X {x}) n ^ ^ 0} > c Cap,__^ (A). (45) 

Proof. The three statements follow from Theorem 2.1 and Proposition 4.3 
appUed respectively to the stochastic process U, U{t) = {u{t,x),x G M^} 
with t e I, and U{x) = {u{t,x),t e [to,T]} with x e J. Notice that by (40) 
and (41), the parameters 7 and a in Theorem 2.1 are 7 = , a — 2 — P, 
and m — k + 1, m — k, m — 1, respectively. □ 



Remcirk 4.2 Since the probability of visiting translates of a compact set A 
decreases to as the distance of this translated set to the origin tends to 
infinity, it is not possible to replace [—N,N]'^ by M.'^ in the above theorem. 
In contrast, this will be possible in the upper bounds of the next theorem. 

Theorem 4.5 Assume (H^). Let I, J be compact subsets of [to,T] and 
M*^, respectively, each with positive Lebesgue measure. Then, the following 
properties hold. 

(1) There exists a positive constant c — c{I, J, (3,k,d) such that, for any 
Borel set A C R"^, 

P{u(I X J) n A ^ 0} < c n^_2(k+ix{A). (46) 

2-/3 

(2) For any t & I, there exists a positive constant c — c{J, (3,k,d,t) such 
that, for any Borel set A C W^, 

P{u {{t} X J) n A ^ 0} < c n,_^ (A). (47) 

(3) For any x & J, there exists a positive constant c — c{I, (3,k,d,x) such 
that, for any Borel set A C W^, 

P{u (/ X {x}) n A ^ 0} < c H,__^ (A). (48) 

Proof. We first note that if we replace d in the Hausdorff dimensions of the 
bounds by any 7 e ]0, d[, then these statements would be a consequence 

of Theorem 2.4 applied respectively to the stochastic processes U, U{t) = 
{u(t,x),x e M^} with t e I, and U{x) = {u{t,x),t E [to,T]} with x e J. 
Indeed, assumption (1) of Theorem 2.4 is given in (42). Moreover, since 
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U is a. Gaussian process, the right-hand side of (33) yields the vahdity of 
hypothesis (2) of Theorem 2.4, with 5 = 

The improvement to 7 = d is obtained by applying Theorem 2.6 to each 
of the stochastic processes mentioned before. Let us argue with the process 
U for the sake of illustration. From (33), we easily deduce that 

< E [exp(cX)] = C, 

where X stands for a standard Normal random variable. Thus, when m — 
k + 1, 6 = ^y^, the left-hand side of (25) is bounded by a constant times 

4(fc+l) 

the square of the volume of Rj, that is, Cs 2-/3 . Hence, the assumptions of 
Theorem 2.6 are satisfied. 

The proof of the theorem is complete. □ 
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